Charge ordering in 0-(BEDT-TTF) 2 RbZn(SCN) 4 : 
Cooperative effects of electron correlations and lattice distortions 
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Combined effects of electron correlations and lattice distortions are investigated on the charge 
ordering in <9-(BEDT-TTF) 2 RbZn(SCN) 4 theoretically in a two-dimensional 3/4-filled extended 
Hubbard model with electron-lattice couplings. It is known that this material undergoes a phase 
transition from a high-symmetry metallic state to a low-symmetry insulating state with a horizontal- 
stripe charge order (CO) by lowering temperature. By means of the exact-diagonalization method, 
we show that electron-phonon interactions are crucial to stabilize the horizontal-stripe CO and to 
realize the low-symmetry crystal structure. 

PACS numbers: 71.45.Lr, 71.10.Fd, 63.20.Kr, 71.30.+h 



I. INTRODUCTION 

Two-dimensional (2D) strongly correlated electron sys- 
tems have attracted much attention in condensed mat- 
ter physics. For bis(ethylenedithio)-tetrathiafulvalene 
(BEDT-TTF) salts, which are 2D strongly cor- 
related organic conductors with various molecular 
arrangements^ many theoretical scientists use the ex- 
tended Hubbard model in order to understand the role 
played by long-range Coulomb interactions. It should 
be noted that not only electron-electron interactions but 
also electron-phonon couplings can be significant because 
the crystal structure is often altered at a phase tran- 
sition. These materials exhibit many interesting phe- 
nomena such as metal-insulator (M-I) transitions with 
charge ordering, the realization of a spin liquid stated 
superconductivity, 4 and so on. In particular, transitions 
to states with a charge order (CO) have been investigated 
theoretically£&£&2 and experimentally! 10 : 11 

#-(BEDT-TTF) 2 RbZn(SCN) 4 (called (9-RbZn for sim- 
plicity hereafter) is a representative compound which 
undergoes a CO transition accompanied with a struc- 
tural deformation by lowering temperature ! ' The 
ground state of the (9-RbZn salt is an insulator with 
the horizontal-stripe CO along t P 4 bonds (HCO-i p4 ) 
shown in Fig. [TJb), where sites 1 and 4 are hole-rich, 
as confirmed by the X-ray structural analysis? 1 ^ 13 C- 
NMR measurement o 13 : 14 and polarized Raman and IR 
spectroscopy! 15 ' 16 The importance of long-range electron- 
electron interactions is well recognized and the mecha- 
nism for stabilizing the HCO-i p4 has been argued mainly 
on the basis of the low-symmetry structure. Because the 
lattice distortions are coupled with the electron system, 
giving rise to the first-order transition sensitive to the 
crystal structure, electron-phonon interactions are also 
important. 

Quite recently, Iwai et al. have observed photoinduced 
melting of CO in the 6>-RbZn salt and in a- (BEDT- 
TTF) 2 13 (called a-l3 for simplicity hereafter) by fem- 
tosecond reflection spectroscopy* 1 ^ The (9-RbZn and a-l^ 



salts show large and small molecular rearrangements, re- 
spectively, at the M-I transition. Their photoinduced dy- 
namics are qualitatively different: the 0-RbZn salt shows 
local melting of CO and ultrafast recovery of CO irre- 
spective of temperature and excitation intensity, while 
the (X-I3 salt shows critical slowing down. Thus, it is 
important to show how electron-phonon interactions are 
significant in the #-RbZn salt. 

This paper is organized as follows. After a brief expla- 
nation of the model in the next section, we present the 
exact-diagonalization results in Sec lIIII for the hole-hole 
correlation functions, the hole densities, and the modu- 
lations of transfer integrals, from which we propose that 
the low-symmetry structure of the #-RbZn salt at low 
temperature is reproduced by introducing electron-lattice 
couplings in the model based on the high-symmetry 
structure at high temperature. Then, in Sec lIVl we dis- 
cuss the numerical results on the basis of a perturbation 
theory from the strong-coupling limit. A brief summary 
is given in SecfVl 



II. MODEL 

We start with the following extended Hubbard model 
at 3/4- filling with electron-phonon couplings of transfer- 
modulation type, 

+ V C 53 n i n 3 + V P 53 nin i '• •) 
(ij)c {i,j) P M=c,pl,p4 



where c\ a creates an electron with spin a at site i, 

n i,a = c i.cr C i,o-j and U, = Tlj.g. tjj = i^F^ , t?^ or 

£p 4 T is the transfer integral for the bond between the z-th 
site and its nearest-neighbor j-th site along the c, pi or 
p4 bond. (i,j) c & n d {i,j)p denote the nearest-neighbor 
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BEDT-TTF salts are estimated from the extended 
Hiickel calculation^ We suppose here from the exper- 
imental data^ that tf 1 = 0.035, t^ T = 0.095, and 
ip 4 T = —0.095 correspond to the high-symmetry struc- 
ture at high temperature, where the pi and pA bonds are 
equivalent. 

We use the exact-diagonalization method for elec- 
trons, regard phonons as classical variables, and ap- 
ply the Hellmann-Fcynman theorem to impose the self- 
consistency condition on the phonons (y c , y p \, y p i) de- 
termined by 



FIG. 1: (Color online) Anisotropic triangular lattice for 6- 
RbZn salt: (a) high-symmetry structure at high tempera- 
ture, and (b) low-symmetry structure at low temperature with 
horizontal-stripe CO along t V 4 bonds. The thin and thick cir- 
cles represent the hole-poor and the hole-rich sites, respec- 
tively. 



pairs i and j along the c bond and the p bond, respec- 
tively. U represents the on-site Coulomb interaction, and 
V c ( p ) is the intersite Coulomb interaction between the i- 
th site and the j-th site on the c(p)-bond. m is the i- 
th molecular translation or rotation explained later from 
the equilibrium position in the high temperature phase. 

and are the electron-phonon coupling strength 
and the elastic coefficient, respectively. For simplicity, 
we perform variable transformations as 



a^u^ = y^ , — ^ = s M . (2) 



III. NUMERICAL RESULTS 

We show energy values in the unit of electron volt (eV) 
in the following. For electron-phonon couplings, we con- 
sider three kinds of molecular displacements: translations 
in the c-direction (u c ), those in the a-direction (it p i), and 
rotations (u P 4) as discussed in terms of elevation angles 
by Watanabe et al^ We suppose that translations in 
the c(a)-direction contribute to modulations of the trans- 
fer integrals on the cl-(pl-) and c2-(p3-)bonds. We also 
assume that molecular rotations produce differences be- 
tween the transfer integral on the p2-bond and that on 
the p4-bond. Then, the modulated transfer integral on 
each bond reads 
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(3) 



The signs here are so determined that y M > corresponds 
to a deviation from the high-temperature crystal struc- 
ture toward the low-temperature one of the #-RbZn salt. 
In this section, we adopt U — 0.7 and V c , V p < 



U/2 
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Typical values for transfer integrals in 



We vary mainly coupling strengths, s c , s p i, and s P 4. 



A. Case without electron-phonon couplings 

At the beginning, we clarify the ground-state prop- 
erties of the model ([1} in the absence of lattice distor- 
tions. In this case, the exact-diagonalization studies of 
the 4x4-site cluster show that the hole densities are uni- 
form {{n'l) = 1 - (n t ) = 0.5) in any combination of V c 
and V p (not shown). Degenerate CO states are mixed 
owing to the finite-size effect. Then, we have calculated 
correlation functions in this cluster. 

Figure [5] shows the hole-hole correlation functions 
without electron-phonon couplings on the V c — V p plane. 
( n i n i+i(i+2)) m denotes the hole-hole correlation function 
along the /i-bond between the i-th and its neighbor (its 
second-neighbor) site, {n^n 1 !-) M being nearly equal to 0.5 
means that ({n^), (rij)) = (1 — S, 1 — S) is equally mixed 
with ((5, S) with small 5. On the other hand, very small 
(n^n^) M represents that ((n^ 1 ), (n 1 -}) = {1 — 5, S) is mixed 
with (6, 1 — 6). Namely, (n^rij) is close to (0.5) when 
the charge disproportionation between the z-th and the 
j-th sites is large (small). 

For large V c , e.g., for (V c , V p )=(0.35, 0.01), 
(«!> c ~ 0.05, (n*r$ +2 ) e ~ 0.44, {n^ +1 ) p ~ 0.23, 
and {n r ln^ + 2) p ~ 0.39. This result means that the 
diagonal-stripe CO is the largest at this point among 
all hole-hole correlations. On the other hand, for small 
V c , e.g., for {Vc, V P )=(QM, 0.35), (n^ +1 ) c ~ 0.46, 
(»?»i+a>c ~ 0.46, (n?nj +1 ) p ~ 0.03, and (n^ +2 ) p ~ 
0.46. Here, the vertical-stripe CO is the largest correla- 
tion. In the ground state with comparable V c and V p , the 
diagonal- and the vertical-stripe COs coexist. Regarding 
the steep variation in the contour plots of the correlation 
functions in Fig. [2] as a phase boundary, we obtain the 
ground-state phase diagram on the V c -V p plane in Fig. [3] 
The phase diagram consists of three uniform phases with 
different hole- hole correlations; the coexistent phase (I), 
the phase with dominant vertical-stripe CO correlation 
(I), and the phase with dominant diagonal-stripe CO 
correlation (IE). This phase diagram is consistent with 
the previous worki^ The most important thing here is 
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FIG. 2: Hole-hole correlation functions as a function of V c and V p for U = 0.7 without electron-phonon couplings: (a) (nj 1 nj +1 ) c 
and (b) (n^n^ +2 ) c are correlation functions along vertical bonds, (c) (nfnj +1 ) p and (d) (n.^nj +2 ) p are correlation functions 
along diagonal bonds. 



that we cannot find the horizontal-stripe CO phase in 
any combination of V c and V p , which is observed experi- 
mentally in the (9-RbZn salt at low temperature. 

Below we mainly use V c = 0.31, V p = 0.27 in phase I 
of Fig. [3l where this ratio of V c /V p is regarded as ap- 
propriate for the 0-RbZn salti^ To check whether these 
values themselves are appropriate, we have calculated the 
optical conductivity spectra by means of the continued 
fraction expansion to show them in Fig. |4p^ It is found 
that for both polarizations E they have a broad peak 
around w~0.4 eV and the conductivity for E\\a is larger 
than that for E\\c, which are consistent with the experi- 
mental findings. 25 Therefore, these interaction strengths 
are quite reasonable. 



B. Case with electron-phonon couplings 

Each of the three kinds of electron-phonon couplings is 
studied at the particular point of V c = 0.31 and V p = 0.27 
by the exact diagonalization. Because we need to ob- 
tain the displacements consistently with bond densities, 
we used the smaller system of the 12 sites hereafter and 



compared some results with those of the 16-site cluster. 
The conclusion is found to be unchanged and, in fact, un- 
derstood from the perturbation theory from the strong- 
coupling limit, as discussed later. 

First, we see the effect of molecular rotations (sp4). 
The s p 4-dependence of the modulations of transfer in- 
tegrals and the hole densities is shown in Fig. E^a) 
with s c = s p i = fixed. For small s P 4, as all the 
phonons are undistorted and the hole densities are uni- 
form ((n^) = 0.5), the ground state remains the coex- 
istent state of the vertical-stripe CO and the diagonal- 
stripe CO. With increasing s P 4, the system discontin- 



uously changes at a critical point s 



pi 



0.04, to the 



CO state of broken symmetry. This CO pattern is the 
HCO-i p4 ((n£), {n'D > (n§), (n^)), which agrees with 
the experimental findings^ Therefore, the coupling with 
molecular rotations (s p 4,) plays an important role to real- 
ize this HCO state. It is noted that y P 4 increases almost 
linearly after the critical point where the HCO-i P 4 state 
is stable. This discontinuous change at the critical point 
would be caused by the finite-size effect: the undistorted 
state is stabilized by quantum tunneling between differ- 
ent CO states, but its energy gain would vanish in the 
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FIG. 3: Phase diagram for U = 0.7 without electron-phonon 
couplings. All phases are uniform, but the dominant hole- 
hole correlations are different. The vertical-stripe CO and 
the diagonal-stripe CO coexist in phase I. The vertical-stripe 
CO is the largest correlation in phase H, while the diagonal- 
stripe CO is the largest correlation in phase HI. 



— E//a 




oj (eV) 

FIG. 4: (Color online) Optical conductivity spectra of the 
model {1} on the 4x4-site cluster for electric fields, E\\a (thin 
line) and E\\c (thick line). We use U = 0.7, V c = 0.31 and 
V p — 0.27. The 5-functions appearing in the continued frac- 
tion expansion are broadened with width 77=0.01. 



thermodynamic limit. In fact, y P 4 rises up linearly by 
infinitesimal s p 4 in mean-field calculations^ so we guess 
that the critical point approaches zero in the thermody- 
namic limit. 

Next, we consider the effect of translations in the 
c-direction (s c ). For large s c , they also stabilize the 
HCO-t p 4 state as shown in Fig. G3h). In the case of 
s c , the charge disproportionation between the hole-rich 
sites ((nj - ), (nj)) and the hole-poor sites ((n^), (^3)) is 
smaller than the case of s P 4, but the translations in the 
c-direction further stabilize the HCO-i^. 

The situation regarding s p i is quite different from the 
cases of s p a and s c . The s pi -dependence of the transfer 
modulations and the hole densities at four sites is shown 
in Fig. [He). In the presence of only y p \ , which represents 



translations in the a-direction, the CO pattern is still 
a horizontal-type but the holes are localized on the t p \ 
and t; P 3 bonds. Thus, the experimental findings are not 
reproduced if only s p \ exists. 

From Figs. [5ja)-[5fc), we conclude that the effects of 
molecular translations in the c-direction (s c ) and molec- 
ular rotations (s p /C) are stronger than that of translations 
in the a-direction (s p i) in this salt. We also come to the 
same conclusion in the other phases of Fig. [3J In Fig. 
IH1 we perform the same calculations (except for the pa- 
rameters shown) as in Fig. [5ja) for phases II and II. The 
s P 4-induced HCO-7j p 4 state is stable in both of phases II 
and III when s P 4 is large enough. Therefore, this char- 
acteristic phonon-induced feature is widely obtained on 
the V c -V p plane, including V c = V p = (irrespective of 
hole-hole correlations in the uniform ground state with- 
out electron-phonon couplings). 

Finally, we simultaneously take the three kinds of dis- 
placements into account. Considering a situation with 
thermal- or pressure-induced expansion or compression, 
we vary intersite Coulomb interactions. In Fig. [3 we 
show the modulations of transfer integrals and the hole 
densities as a function of nearest-neighbor Coulomb inter- 
actions. For small r, only y p \ is finite with y c and y p 4 be- 
ing zero. This situation is quite similar to Fig. [5tc), and 
hence the most stable state is the HCO-r;pi&;7; p 3. As r ap- 
proaches unity, which is a realistic value for 0-RbZn, the 
experimentally observed HCO-i P 4 state becomes more 
stable than the HCO-i p i&i P 3 by gradually increasing y c 
and y P 4 and rapidly decreasing y p \. In the vicinity of 
r = 1 for s c = 0.05, s p i = 0.17, and s P 4 = 0.06, y c 
and y P 4 approximately correspond to the low-symmetry 
structure of #-RbZn at low temperature. 

From the X-ray structure analysis ji 2 - the molecular 
translations in the c- and a-directions are found to be 
u c ~ 0.13 A and u p x ~ 0.15A, and the corresponding dis- 
tortions are y c = 0.0185eV and y pl = 0.0735eV. They 
give a c - 0.14eV/A and a pl ~ 0.47eV/A. The present 
coupling strengths s c = 0.05eV and s p i = 0.17eV lead 
to K c - 0.38eV/A 2 and K pl ~ 1.32eV/A 2 . A rough es- 
timation of phonon frequencies by oj^ = -J K^/m with 
m being the reduced mass of two BEDT-TTF molecules 
shows lo c ~ 23cm" 1 and u> p i ~ 43cm~ 1 , which are quite 
reasonable values.— 

As r further increases beyond two, all modulations 
approach zero. In particular, y p \ decreases monotoni- 
cally. We explain these behaviors with a perturbation 
theory from the strong-coupling limit in the following 
section. Thus, electron-phonon interactions are crucial 
to stabilize the horizontal-stripe CO and to realize the 
low-symmetry structure of 0-RbZn at low temperature. 
All these results are consistent with mean-field studies of 
the same model. 26 
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FIG. 5: (Color online) (a) s P 4-dependence of modulations of transfer integrals (left) and hole densities (right) for s c =s p i=0, 
(b) their s c -dependence for s p i=s P 4=0, and (c) their s p i-dependence for s c =s P 4=0, with U = 0.7, V c = 0.31 and V p = 0.27. 



IV. DISCUSSION 



site read 



In this section, we consider why the HCO-i p 4 is the 
most stabilized by the electron-phonon couplings. It is 
easily understood in the strong-coupling limit. We as- 
sume that, in this limit, the holes are perfectly localized 
on the t P 4 bonds, namely (nf) = (roj) = 1 as shown 
in Fig. Hlb). In the ordinary perturbation theory, the 
second- and third-order contributions to the energy per 
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(6) 



where Si is the 5=1/2 spin operator at the i-th site 
on the t p 4 bond. In the strong-coupling limit, the one- 
dimensional half-filled chain is formed along the t pi bonds 
and (Sj • S»+i) = -ln2 + 1/4 ~ -0.443 is the exact 
ground-state energy of the isotropic 5=1/2 Heisenberg 
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FIG. 6: (Color online) s P 4-dependence of modulations of transfer integrals (left) and hole densities (right), (a) in phase H 
(V c = 0.01, V p = 0.35), and (b) in phase JSL (V c = 0.35, V p = 0.01), for U=0.7 with s c =s p i=0 fixed. 



chain. According to the third term in Eq. there is 
energy gain from spin fluctuations, so that the ground- 
state energy of the HCO-i P 4 state is further lowered by 
increasing \t p i\. On the other hand, because all terms in 
Eq. © are positive, E 3 represents energy loss. When \tp±\ 
becomes large, |t c i| becomes small to reduce the energy 
loss. From the relations ([3|), |t C 2| becomes large and \t P 2\ 
becomes small. 

Equations ([5]) and ^ can be rewritten as 

E% oc a + a i p4 T y P i , (7) 
E 3 oc b + c 2 tf T y p4 + c 3 t^J y c , (8) 

where a, b and c$ (i = 1,2,3) are independent of y lll 
and tp 4 T < < t^ T . Note that we ignore in deriving 
the above equations, and the coefficients Cj are positive 
under the condition 2V C > V p . We can infer from Eq. ([7]) 
that y P 4 increases linearly as a function of s P 4 when the 
ground state is the HCO-t P 4 state. Because of the high 
symmetry with respect to y p \ <-> —y p i, Eqs. |(7J| and 
([8|) are independent of y p \ . Once the distortions y P A and 
y c are substantially large, however, the distortion y p \ is 
numerically found to stabilize the HCO-£ P 4 further. 

Before closing this section, we make a brief comment 
on the spinless fermion case^i In the U — > oo limit of 
the model([T]), the HCO state becomes unstable even if 
the modulations of the transfer integrals are introduced 



because the energy gain from spin fluctuations is absent. 
Thus we demonstrate that cooperative effects of electron 
correlations and electron-phonon couplings are important 
in the 0-RbZn salt. The stability of the HCO state rel- 
ative to the diagonal CO state suggested in the spinless 
fermion case without electron-phonon couplings 27 might 
be caused by the absence of phase factors, which are 
present if the U —> oo limit is naively taken. 



V. SUMMARY 

We have investigated the cooperative effects of elec- 
tron correlations and lattice distortions on the charge 
ordering in the #-RbZn salt. By means of the exact- 
diagonalization method for systems with up to 16 sites, 
we have calculated the hole-hole correlation functions, 
the hole densities and the modulations of transfer inte- 
grals to clarify the role of electron-phonon couplings in 
this salt. 

In the absence of electron-phonon interactions, there 
appear three uniform phases characterized by quantum 
tunneling between the vertical-stripe COs, that between 
the diagonal-stripe COs, and that between the vertical- 
and the diagonal-stripe COs. It is found that all of 
these phases are changed into the horizontal-stripe CO 
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by the introduction of electron-phonon couplings relevant 
o Yc to the 6*-RbZn salt. We can mostly reproduce the low- 

A Yp1 symmetry structure of 0-RbZn at low temperature by us- 

□ Yp4 ing the model based on the high-symmetry structure at 

high temperature and by choosing the coupling strengths 
appropriately. 

We conclude that the structural deformation assists 
the horizontal-stripe CO that is experimentally founds 
In particular, the effects of molecular translations in 
the c-direction and molecular rotations are found to be 
stronger than that of translations in the a-direction. 
Thus, electron-phonon couplings are significant to stabi- 
lize the HCO-i P 4 and to realize the low-symmetry struc- 
ture of 6*-RbZn. Otherwise the long-range Coulomb inter- 
actions favored a different CO pattern. With the help of 
the perturbation theory from the strong-coupling limit, 
we easily understand the mechanism for stabilizing the 
HCO-t P 4 state and the linearly increasing y p 4 and y c 
modulations. 



FIG. 7: (Color online) Three kinds of transfer modulations 
and hole densities at four sites for U = 0.7, V c = 0.31r 
and V p = 0.27r. Coupling strengths are fixed at s c = 0.05, 
s p i — 0.17 and s P 4 = 0.06. The solid (dashed-dotted) line 
in the upper panel represents the magnitude of the molec- 
ular translations in the c(a)-direction in the low-symmetry 
structure of #-RbZn at low temperature.— The dashed line 
represents the magnitude of the molecular rotations in this 
structure. 
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